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Abstract 

A diffusion molecular channel is a channel with memory, as molecules released into the medium 
hit the receptors after a random delay. Coding over the diffusion channel is performed by choosing 
the type, intensity, or the released time of molecules diffused in the environment over time. To avoid 
intersymbol interference (ISI), molecules of the same type should be released at time instances that are 
sufficiently far apart. This ensures that molecules of a previous transmission are faded in the environment, 
before molecules of the same type are reused for signaling. In this paper, we consider ISI-free time- 
slotted modulation schemes. The maximum reliable transmission rate for these modulations is given 
by the constrained coding capacity of the graph that represents the permissible transmission sequences. 
However, achieving the constrained coding capacity requires long blocklengths and delays at the decoder, 
making it impractical for simple nanomachines. The main contribution of this paper is to consider 
modulations with small delay (short blocklength) and show that they get very close to constrained 
coding capacity. 
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I. Introduction 


Molecular communications (MC) via diffusion is one of the most promising approaches for 
communications among nanonetworks [[T|, Q. In this approach, information is conveyed to a 
nano-receiver by nano-transmitter’s choice of the concentration, type, or the release time of 
the molecules diffused into the medium. The motion of the released molecules are described 
by a Brownian motion process, that can be with or without drift. A molecule released from a 
transmitter can follow different trajectories before reaching a receiver. Thus, diffusion-based 
communication suffers from intersymbol interference (ISI) due to molecules from previous 
transmission that follow longer trajectories before hitting the receiver. 

The effect of ISI in diffusion based MC has been studied extensively in the literature Q- 
Q, with the general conclusion that ISI reduces the performance in a communication setup 
consisting of one transmitter and one receiver. The authors in Q, Q have proposed two 
modulation schemes. Concentration Shift Keying (CSK) and Molecular Shift Keying (MOSK), 
where both suffer from the ISI caused by molecules from previous transmissions. Authors in [|^ 
have proposed the Molecular Concentration Shift Keying (MCSK) modulation scheme, where the 
main idea is to use the distinct molecule types for consecutive time slots at the transmitter, thus 


effectively suppressing the ISI. In [10| the authors take advantage of a limited memory in the 
transmitter to propose an adaptive transmission rate scheme that reduces ISI. An on-off MOSK 
modulation scheme has been proposed in pT| which increases the transmission rate compared 
to MOSK but still suffers from the ISI. In p^ , a modulation, called run-length hybrid aware, 
has been presented, which considers the runs (the same value occurs in several consecutive 
bits) for encoding: the run-value is encoded in molecule type and the run-length is encoded by 
CSK. Though suffering from the ISI, this scheme is stated to improve the transmission rate in 
comparison with CSK and MOSK considering the ideal channel with no transmission errors. 


Authors in [ 131 propose a pre-equalization method to mitigate ISI where the signal at the receiver 


DRAFT 


August 12, 2015 




3 


is considered as the differenee between the number of the reeeived molecules of eaeh type while 
two types of moleeules is used at the transmitter. Further, several ehannel eoding sehemes are 


proposed to mitigate ISI in diffusion based MC [14|-[16| where they employ additional bits to 
intended input eodewords in order to eompensate errors from the ISI or noisy ehannel. 

Motivation for ISI-free modulations: Given the generally negative effeet of ISI, it is of interest 
to design modulation or ehannel eoding sehemes that would prevent or mitigate the ISI. To 
mitigate the ISI, the transmitter needs to make sure that moleeules of the same type are released 
at time instances that are suffieiently far apart. Therefore, the transmissions will not have any 
superposition. This simplifies the reeeiver, who monitors the presenee and intensity of various 
moleeule types in the environment. Furthermore, ISI-free modulation sehemes are more robust 
when the reeeiver knows little about the ehannel state information (CSI). A high reeeiver moleeule 
intensity ean be either due to a high transmission amplitude or a very eonductive ehannel. In 
laek of CSI and potential presenee of ISI, observing a eertain intensity of a moleeule type, the 
reeeiver has to meet the ehallenge of deeiding whether the reeeived moleeule density is due to 
superposition of various suecessive transmissions or due to temporary ehannel eonductivness. 

In this paper, we study ISI-free modulation sehemes for time-slotted transmission model. The 
problem of finding the maximum aehievable rate is eonverted into the design of a eodebook, 
eonsisting of string of symbols, with eaeh symobl indieating the type and eoneentration of a 
moleeule transmission. The ISI-free eondition implies that symbols of the same type should not 
be used elose to eaeh other. This problem ean be expressed in terms of the eonstrainted eoding 
problem. Constrained eoding has found applieations in storage and eommunieation systems; its 
basic idea is to avoid transmission of input sequenees that are more prone to error (e.g. see CD 
Chapter 1]). Similarly, in our setup, we would like to avoid input sequenees that eause ISI. 

Aeheiving the eonstrained eoding eapaeity requires bloeklengths tending to infinity and eom- 
plieated eoding sehemes. On the other hand, due to the limited resourees available at nano-level, 
only simple sehemes are eligible for being implemented in praetiee. Further, we desire short 
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blocklengths (limited delay) and modest memory resources. Our main contribution in this paper 
is to construct simple ISI-free modulation schemes for a given maximum transmission delay. It 
is shown (numerically) these simple schemes are near optimal, i.e., they become close to the 
constrained coding capacity. 

The rest of this paper is organized as follows. In Section ISI-free modulation is formally 
defined and modeled as a constraint graph. A summary of the notation used in this paper is given 
in Table 1^ at the beginning of this section. Section [ni| introduces the ISI-free modulation capacity 
under the limited delay (blocklength). In Section a class of modulation strategies based on 
a modified version of the constraint graph is proposed as a lower bound on the capacity. The 
numerical results and conclusions are provided in Section |V] and |V^ respectively. 

II. ISI-FREE MODULATION PROBLEM FOR DIFFUSION CHANNELS 

Herein, we study ISI-free modulation schemes for the following setup: time is divided into 
slots of a fixed duration. The transmitter operates at the beginning of each time-slot: it either does 
not release any molecules, or chooses one molecule type and releases it into the environment. 
For simplicity of exposition, we assume that the information is coded into the molecule type 
only, and not in the number of released molecules (on-off keying modulation); extensions to 
amplitude modulation is also possible. Therefore, the transmission sequence can be expressed 
as a string of symbols, each of which is either a — (no transmission), or Mi when molecule of 
type i is used. 

To avoid ISI at the receiver, there should be some time delay between two consecutive use of 
molecules of the same type. The amount of this delay depends on the communication medium. 
The diffusion channel is assumed to have k time slot memory. Thus, we restrict the transmission 
sequence by requiring that there should be at least k symbols between any two consecutive 


'please note that we use the words “code” and “modulation” and their derivatives interchangeably. 
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TABLE I: Summary of the notation 


Notation 

Description 

k 

Diffusion channel memory in terms of time slot 

N 

The number of molecule types 

Mi 

Symbol corresponding to molecule type i 

- 

Symbol corresponding to no molecule transmission 

Si 

i-th symbol in a symbol sequence 

s 

a symbol string or a state 

1 

the length of a binary string 

m 

length of a symbol string 

d 

Delay (blocklength) constraint 

£ 

Encoding function 

V 

Decoding function 

C{d) 

ISI-free capacity given blocklength d 

TZtid) 

Lower bound rate given blocklength d 

W(s) 

Set of all possible uniquely decodable subsets 


corresponding to edges that exit from state s 

U{s) 

a member of set W(s) 

Uu{s) 

a full prefix free binary code mapped to sequences in U (s) 


symbols M* for some parameter k. For instance, if k = 2, Mi — M 2 Mi and Mi-Mi are valid 

sequences, but Mi — Mi M 2 is invalid. 

Definition 1. A sequence of symbols (si, S 2 , S 3 ,, s„) is called ISI-free if no the symbol of the 
form Mr appears more than once in any consecutive substring of length k + 1. The symbol, 
is however, allowed to be repeated more than once. In other words, for any 1 < i < n — k, the 
substring Si, Sj+i,..., Si+k, has the property that Sj = Sr for j f r ^ {i,i + 1,... + k} only 

if Sj Sr . 
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Mi/1 



Fig. 1 : State diagram of the MCSK modulation scheme 


We assume that the reeeiver has reeeptors that ean deteet the presenee of any moleeule type 
over time. Having avoided the ISI, the reeeiver ean fully reeover the transmission sequenee 
of symbols. Then, the task of the transmitter is to eneode information bits into a transmission 
sequenee that obeys the restrietion given above. 

Example: An example of an ISTfree modulation is the MCSK modulation seheme, proposed 
in Q, for diffusion ehannels with one symbol memory {k = 1)0 MCSK uses the symbol set 
M = {—, Ml, M 2 } and prevents ISI by avoiding the transmission of the same moleeule type in 
two eonsequent time slots. In binary MCSK, ‘0’ input bit is mapped to ’ symbol; to transmit bit 
‘1’, in even time slots symbol Mi is used while in the odd time slots, symbol M 2 is used. Note 
that quaternary and higher modulation orders ean be developed by using more eoneentration 
levels of different moleeule types. Fig. demonstrates the state diagram for the MCSK. In 
this figure, eaeh edge is labeled by x/y, where x and y indieate the transmitted moleeule and 
eorresponding input bit, respeetively. Observe that the MCSK seheme eommunieates at the rate 
of 1 bit per symbol. 

Observe that after transmission of the symbol ‘ — i.e., no moleeule transmission, in faet all 


^ MCSK can be easily extended to the case of general channel with k memory by using A: + 1 different molecule types. 
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Fig. 2: State diagram of ISI-free sequences of the symbol set M = {—for channel with k = 1 


of the three symbols Mi and M 2 and ‘ —’ can be transmitted without arising ISI. This degree 
of freedom is depicted in the state diagram of Fig. If we move on the edges of this graph 
and write down the sequence of symbols written on the edges, we get an ISI-free transmission 
sequence. Conversely, any ISI-free sequence corresponds to a path on this graph (we allow an 
edge to occur arbitrarily many times in a path). Therefore, the number of paths of length m 
on this graph give us with the maximum number of messages that we can send by transmitting 
m symbols. A transmission rate R is equivalent with sending a binary sequence of length mR 
using m symbols, i.e., the number of messages is 2"*^ for a code of length m. Therefore, we 
are interested in the logarithm of the number of paths in the graph, divided by m. 

When the channel memory k and the number of molecule types are arbitrary, one can draw 
a state diagram similar to Fig. as follows: 

Definition 2 (Constraint graph). Each vertex (state) is defined as the vector of the molecule 
types of the last k transmitted symbols and each of the edges emanating from a vertex (state) 
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represents an allowed symbol for transmission given the last k transmitted symbols. In other 
words, the nodes of the graph are k-tuples (si, § 2 , • • • , Sk) where Si = Sj for i j occurs only 
if Si = Sj = There is a directed edge from vertex (si, S 2 , • • • , su) to (s 2 , S 3 , • • • , s^, s^+i), 
only if {si, S 2 , ■ ■ ■ ,Sk,Sk+i) is an ISI-free sequence, i.e., either Sk+i =‘—’, or Sk+i f si. The 
edge from vertex (si, S 2 , • • • , s^) to (s 2 , S 3 , • • • , s*,, s^+i) is labled by the symbol Sfc+i. 


The number of paths of length m on a direeted graph (starting from a given node) is denoted 
by N{m) and is a simple eombinatorieal problem that ean be found using reeursive equations. 
The logarithm of the number of paths of length m divided by m, as m eonverges to infinity, 
ean be readily found in terms of the adjaeeney matrix of the graph, as it is a speeial ease of the 
classical problem of constrained coding Theorem 3.4]. 


Theorem 1. ( [17 Theorem 3.4]) Let Q be an irreducible constraint graph representing an 


irreducible constrained system. Then, 


lim sup — log A^(m) = log 2 \{Ag), 

n—>-oo m 


where Ag is the adjacency matrix of the graph Q and X{Ag) is the largest of the absolute values 
of eigenvalues of Ag. 


The connection between constrained coding and ISI-free coding is known and studied in 
applications such as magnetic recording constrained coding [[^ Chapter 1]. The earliest results 


on constrained coding problem are due to Shannon [1^ who considered teletype and telegraphy 
systems as two examples of channels with restrictions on the input symbol sequences. A general 
result for transmission of a binary symmetric source over a finite state discrete noiseless channel 
is given by Shannon in p^ . 

According to Theorem one can easily verify that the maximum achievable transmission rate 


DRAFT 


August 12, 2015 




9 


for the graph of Fig. is C = 1.2716. The adjaeeney matrix for this graph would be as follows: 

- Ml M2 

Ml 
M2 

and for the given adjaeeney matrix \{Ag) = 2.4142, so the maximum aehievable rate equals 
log2(2.4142) = 1.2716. But aehieving this rate requires long bloeklength that implies long 
memory and transmission delay at the transmitter and reeeiver. Nonetheless, we now show that 
one ean reaeh the transmission rate of 1.25 with bloeklength (delay) of one. We do this by 
slightly modifying the MCSK seheme (MCSK was transmitting at rate 1 bit per symbol). 

Consider the modulation seheme given in Fig.j^with the same x/y eonvention as before. Here, 
for instanee if we are in state —, we send symbol — if the next message bit we want to transmit 
is zero; we send Mi if the next two bits are 10 and we send M 2 if the next two bits are 11. In 
this modulation seheme, the average transmission rate in states —, Mi, M 2 are 1.5,1,1 bits per 
symbol respeetively. Also, the probability of being in states —, Mi, M 2 ean be eomputed using 
eorresponding linear equations over the states, whieh will be equal to 0.5, 0.25, 0.25, respeetively. 
Therefore, this seheme aehieves the rate of 1.25. 

The implieation of the above eonstruetion is the existenee of near optimal modulation strategies 
that are effieient in terms of bloeklength and deeoding delay, as the rate 1.25 is very elose to 
the eonstrained eoding limit of 1.2716.... The main result of this paper is eonstrueting simple 
memory-limited ISI-free modulation sehemes for a given maximum transmission delay that 
aehieve rate elose to the eonstrained eoding limit. 

III. ISI-FREE MODULATIONS WITH LIMITED DELAY AND DECODER BUFFER SIZE 

In the previous seetion, we observed that the ISI-free eapaeity, when there is no restrietion 
over delay, is equivalent to the eonstrained eoding eapaeity. In praetieal modulation sehemes for 
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Fig. 3: A modulation strategy for the diffusion channel and the symbol set M = {—, Mi , M 2 } for channel with 
k=l. 

MC, due to limitations on memory resourees and the faet that large deeoding delays may not be 
allowed in some applieations, we are interested in assigning the input bits to ISI-free sequenees 
of symbols with a limited length, i.e., using eodes with a short bloeklength. 

Definition 3 (Delay (bloeklength) eonstraint). Delay constraint d means that the decoding of 
the i-th bit of the input binary string, bi, can be delayed up to d channel uses. More precisely, 
in using a variable length code for transmission of the input data bits, no codeword of length 
more than d should be used. This implies that upon the transmission of the bit bi, it will be 
decoded within d time-slots. 

An ISI-free eode of length n, with a given bloeklength (delay) eonstraints eonsists of 

• An eneoder funetion, 8, mapping the input bit string { 61 , ■ ■ ■ , to an ISI-free sequenee 
of the symbols (si, • • • ,Sm)- The number of output symbols m may depend on the input 
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bit string. The rate of the code is R = n/E[m] (bits per symbol). 

The maximum blocklength d constraint implies the following structure on the encoding 
function S = { 81 , 82 ,- ■ ■): assuming that bits {hi,b 2 ,--- ,bi^) are already encoded into 
symbol string (si, S 2 , • • • , 8 ^), the encoder 8i takes in new bits ..., 6^-^^ and maps it 
to output symbols (sr^+i, ^^^+ 2 , • • • , -Sri+i) of length less than d, i.e., Tj+i —Ti<d. Observe 
that li is the number of input bits that are encoded in the first i-th steps and r* is the number 
symbols that are transmitted during this step. 

• A decoder function Vi at step i that reverses the encoding operation by taking in symbols 
(si,S 2 , • • • ,-Sri+i) and producing bits where k is the number of input bits 

that are decoded in the i-th step. 

Given d, the ISI-free capacity of the system is defined as the maximum rate that is asymp¬ 
totically achievable: 


C{d) 


n 


max 


n,£i,V^ 


( 1 ) 


It is clear that C{d) < C{d') iid < d'. Further, the maximum achievable rate of the system when 
we do not have any delay constraint, C(cx)), is equal to the constrained coding capacity. Therefore, 
constrained coding capacity is always an upper bound for C{d). In this paper, we propose a 
suboptimal approach to determine the encoder and decoder function for ISI-free modulation 
problem under delay constraints. The associated rate gives a lower bound on C{d), which we 
denote by lZ\_{d). Thus, 

R^d) < C{d) < C{oo). 


Simulation results show that lZ\_{d) is close to C(cx)). for small values of d. As a result, C{d) 
is not far from C(cx)). and the proposed scheme is near optimal. To propose the lower bound, we 
restrict to a class of modulation schemes constructed based on state diagram of ISI-free symbol 
sequences. 
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Remark 1. Although our aim is to study diffusion molecular channels, but the delay-limited 
codes that we construct in the next section can be used in other similar constrained coding 
applications. 


IV. Lower bound on C{d) 

Consider a diffusion channel with k symbol memory, where the symbols are selected from 
the set M = , Mi,M at}, in which N is the number of molecule types available to the 

transmitter. In an ISI-free modulation, no molecule type appears repeatedly, in each sequence of 
output symbols of length fc + iQ As discussed earlier, the ISI-free sequences can be represented 
as a path on the constraint graph of Definition One example of such a graph is given in Figure 

To construct a good code, we used the graph given in Figure reaching the performance of 
1.25 which is close to the optimal 1.2716. 

The main idea of this section is to provide a systematic algorithm for constructing codes 
similar to the one given in Figure from the constraint graph of Figure To do this, we 
first define a modified constraint graph of maximum depth d to represent a practical ISI-free 
modulation family under the memory and delay constraints. We then discuss assigning input bits 
to edges of this modified constraint graph. 

In the modified constraint graph, the vertices are not changed but a permissible sequence of 
maximum length d (instead of a single symbol) is assigned to each edge. 

Definition 4 (Constraint graph of depth d). Let d be a natural number. The nodes of the 
constraint graph of depth d are k-tuples (si,S 2 , • • • ,Sk) of permissible sequences, i.e., Sj = Sj 
for i f j only if Si = Sj = ’. Given any natural number t < d and any ISI-free string 


^The number of molecule types N may be less than k, since we are allowed to use — repeatedly in constructing the ISI-free 
sequences. 
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(si, S 2 , • • • , Sk+t), we draw a directed edge from vertex (si, S 2 , • • • , Sfc) to {st+i, St+ 2 , ■■■ , St+k)- 

This edge is labeled by the string (s^+i, < 5 ^+ 2 , • ■ ■ , -St+fc). 

Different ISI-free modulations are obtained by assigning sequences of binary bits to a set 
of edges emanating from each state. The modulator produces the ISI-free sequence by moving 
over the introduced graph based on the input bit sequence and transmitting the corresponding 
symbols. The demodulator reconstructs the bit sequence by moving on the same graph based on 
the received symbol sequences. In this modulation strategy, the demodulator needs to buffer at 
most the last d + k received symbols, imposing a limited symbol buffer of length d + k at the 
receiver. 

More accurately, given a state s = (si,S 2 ,--- ,Sk), the emanating edges from the state are 
labeled with a string of length t < d, Sk+i, Sk+ 2 , • • •, Sk+t such that (si, S 2 , • ■ ■ , Sk+t is ISI-free. 

We would like to select a subset of these edges and assign binary bits to them. This assignment 
specifies the edge we would choose for moving, based on the input bit string. The only restriction 
is that the set of sequences Sk+i, Sk+ 2 , ■ ■ ■, Sk+t of the edges that we choose must be prefix-free 
(uniquely decodable). This would enable the possibility of recovering the input bit string from 
the set of output strings. 

Example 1. For the symbol set M = { — , Mi, M 2 }, d = 2 and k = 1, the constraint graph of 
depth d has three states —,Mi and M 2 . Given the state s = Mi, the allowed edges emanating 

from this state are labeled by { —, M 2 ,-, —Mi, —M 2 , M 2 —, M 2 M 1 }. From these seven edges 

that leave state Mi, there are two edges (labeled by M 2 and —M 2 ) that go from state Mi to 
state M 2 ; there are two edges (labeled by —Mi and M 2 M 1 ) that go from state Mi to itself; and 

there are three edges (labeled by —, -, and M 2 —) that go from state Mi to —. Among these 

seven edges, the five edges {-, —Mi, —M 2 , M 2 Mi, M 2 —} form a possible uniquely decodable 

sequence set (no members of this set is a prefix of the another). 
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Given a state s, let W(s) denote the set of all possible uniquely deeodable subsets eorresponding 
to edges that exit from state s. To eome up with an explieit modulation, for eaeh state s, the 
eneoder parses the input data bits by a full prefix-free eode and maps eaeh eodeword to 

a sequenee from the f/(s) G W(s). 

Example 2. Consider Example^ the state s = Mi and U (s) = {-, —Mi, —M 2 , M 2 M 1 , M 2 —}. 

Given that U (s) has five sequences, to define for this state, we may pick any Huffman code 
on five symbols, e.g. {0,10,110,1110,1111}, and assign it to the sequences ofU{s) respectively. 

Thus, if we are at state s = Mi and the next input bit is 0, we move along edge -; if the next 

two input bits are 10, we move along —Mi, etc. 

Assuming that U{s) = {si, S 2 , • • ■ , s,.} and we assign a binary Huffman eode with string of 
length ii to sequenee Sj, we must have X]i=i = 1- Furthermore, the probability that we 
take the edge along the sequenee Sj will be as the input is assumed to be a uniform and 
independent string of bits. If the length of Sj is symbols, the average transmission rate (bit 
per symbol) given the state s will be equal to 



( 2 ) 


in whieh Es(£) and Es(m) denote the average number of transmitted bits and symbols given the 
state s. The overall aehievable rate by this seheme for large input bit stream ean be eomputed 


by 


^^p(s)Es(m) 


(3) 


where p(s) denotes the stationary probability of the state s and these probabilities are eomputed 
using reeursive linear equations of states probabilities. 

Lower bound on C{d): Any eoding seheme that satisfies the delay eonstraint d would yield 


a lower bound to C{d). Consider a eode eonstrueted based on the eonstraint graph of depth d; 


this eode will satisfy the delay eonstraint of d beeause all the eodewords on edges of the graph 
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are of length at most d symbols. The optimal strategy is determined by the uniquely deeodable 
sequenee set and full prefix eode given eaeh state that aehieves the maximum average rate (for 
large input bit streams) under limited bloeklength eonstraint,/.e., 




max 

for all s 


EsP(s)Es(m)' 


(4) 


The above optimization problem involves maximizing a ratio. The lower bound TZii^d) ean be 
expressed in an alternative linear form, whieh makes it amenable to dynamie programming 
analysis: let 


G{.R) = max Vp(s)E,(£) - R ■ ^p(s)Es(m) 

s s 

= max p(s) (Es(£) — i? ■ Es(m)). (5) 

S 

Then G{R) is a monotonieally strietly deereasing funetion, with G'(O) > 0. Furthermore, G{7l\_{d)) 
0. Therefore, the value of R\_{d) ean be found as the unique positive root of G{R). An effieient 
algorithm for eomputing G{R) would lead to an effieient algorithm for approximating TZ\_{d), 
by employing the biseetion method for eomputing the root of the strietly deereasing funetion 
G{R). 

Dynamic programming (DP) algorithm: The expression for G{R) in equation Q is in a 
form that ean be solved using DP. The state spaee for the DP are the states s (the nodes of the 
graph). Let us interpret ((7(s),/i{/(s)) as the action that we take at state s. The expression for 
Es(£) — R ■ Es(m) depends only on the aetion (U{s), ^ui^)) that we take in state s; it does 
not depend on p(s). Therefore, Es(£) — R ■ Es(m) ean be viewed as the one-step profit, when 
we take the aetion at the state s. And G{R) will be the overall average profit of 

the DP. The state transition is defined in the natural way: upon taking the aetion (f/(s),/i[/(s)) 
at the state s, an outgoing edge eorresponding to members of U{s) is ehosen aeeording to the 
probability distribution given by the length of the binary strings in /i[/(s); we move along this 
edge to get to a new state. 
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The average profit G{R) can be found iteratively using the standard dynamic programming 
algorithm. We now briefly explain an idea to use the structure of our particular DP to simplify 
its implementation: suppose we in an iteration of the DP algorithm. Fix some state s. In each 
iteration of the DP algorithm, we need to maximize (over all actions) the average one-step profit 
plus the profit of the state will end up with. We can do the search over U (s) exhaustively. Then 


the maximizing /it/(s) can be found using the Geometric Huffman Coding approach [ 19|. To see 
this, observe that the average one-step profit plus the profit of the state we reach can be written 
as follows (using the expressions given in equation @): 


( ^ 2-H - ^ 2-^*c(s,) 


( 6 ) 


where c(sj) is the current state cost of the state Sj that we reach from state s with probability 
2~^\ The goal is to maximize @ over all natural numbers {ii, £ 2 , ■ ■ ■, £r) satisfying 2“^* = 1. 
But equation @ can be expressed as the KL divergence —D{p\\q) for p{i) = 2“^* and q{i) = 
2 Rm,-c{si)^ Observe that g(-) is not necessarily a probability distribution. Maximizing —D{p\\q)) 


or equivalently minimizing i9(p||g) over diadic pmf p is the problem considered in [19| and can 
be solved using the Geometric Huffman Coding approach. 


V. Numerical results for the diefusion channel 

In this section numerical results are presented, indicating the near-optimality of the proposed 
modulation scheme. The numerical results illustrates how the system parameters affect the 
achievable bit rates. The parameters investigated here are the number of molecule types exploited 
in the system N, the channel memory k, and the maximum tolerable decoding delay d. 

In Fig. Q the achievable rates are plotted versus the constraint graph depth (i.e. delay d). 
The figure compares the maximum achievable rate for a delay d. The curve is obtained by the 
dynamic programing (DP) algorithm. This figure is plotted for diffusion channels with up to 
three symbols memory, i.e. k = 3. For each channel, the number of molecule types are set to 
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Constraint Graph Depth (Symbols) 


Fig. 4: The effect of the blocklength (constraint graph depth) on the maximum achievable rate. 


k + 1, i.e, = fc + 1. As we see in the figure, small values of d are enough to approaeh the 
eapaeity. 

In Fig. the effeet of number of moleeule types on aehievable rates is studied. Similar to 
Fig. Q the aehievable rates are obtained by the DP algorithm and are plotted for ehannels with 
up to three symbols memory k = 3. For all the plots in this figure the depth of the eonstraint 
graph is set to a single symbol, d = 1. It is observed that as the number of molecule types 
increases, the achievable rate increases. It should be noticed that in Fig. unlike the Fig. 
channels with higher memory have larger rates. The reason is that, the number of molecule types 
are not the same for all the plots as the number of molecule types N increases with k. 

In Fig. the achievable rates are plotted versus the number of molecule types but for different 
values of delays d. As we see in this figure by increasing the depth from one to three, we can 
achieve the capacity for any number of molecule types but the gap occurs starting around N = 5. 
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Fig. 5: Maximum achievable rate vs number of molecule types for d = 1 and fc = 1, 2,3. 


Remember that MCSK was an ISI-free modulation for A; = 1 and N = 2 from the symbol 
set {—, Ml, M 2 }. Table 1^ shows the aehievable rate of different ISI free modulation sehemes 
for k = 1, N = 2 and different values of d. As d inereases the optimal modulation found over 
eonstrained graph gives higher rate very elose to the eapaeity. Observe that the rate obtained for 


d = 2 is the same as that for d = 1. Table III shows the optimal eode that one gets for the ease 
of d = 2, aehieving the rate 1.25. 


VI. Conclusion 

We eonsidered the problem of effieient rate ISI-free modulation in moleeular eommunieation 
under limited delay eonstraints. The ISI-free eondition was expressed by a eonstrained graph 
where maximum aehievable rate of modulations eonstrueted over this graph, i.e, ISI-free eapaeity, 
is well known in eonstrained eoding literature. Motivated by limited resourees in nanomaehines. 
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Fig. 6: Maximum achievable rate vs number of molecule types for fc = 1 and d = 1,3. 


TABLE II: Achievable rate for different ISI fee modulation schemes and different constraint graph depths, k = 1, 
N = 2 and different d values. 


Modulation scheme 

Achieved Rate 

MCSK 

1 

Optimal over constraint graph with d = 1 

(Modified MCSK) 

1.25 

Optimal over constraint graph with d — 2 

(Table III 

1.25 

Optimal over constraint graph with d = 3 

1.2604 

Optimal over constraint graph with d = 4 

1.2617 

Optimal over constraint graph with d = 5 

1.2640 

Capacity for d = oo 

1.2716 
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TABLE III: Optimized modulation for transmitting a uniform binary string over a diffusion channel with one 
symbol memory, according to a constraint graph with depth 2 and the symbol set —, Mi, M 2 ■ 


State: 

- 

State: Mi 

State: M 2 

Symbol Seq. 

Bit String 

Symbol Seq. 

Bit String 

Symbol Seq. 

Bit String 

1 

0,0,0 


0,0 


0,0 

Mi,- 

0,0,1 

M 2 ,- 

0,1 

Mi,- 

0,1 

M 2 ,- 

0,1 

-,Mi 

1,0,0 

-,Mi 

1,0,0 


1,0,0 

-,M2 

1,0,1 

-,M2 

1,0,1 

-,M2 

1,0,1 

M 2 , Ml 

1,1 

Ml , M 2 

1,1 

Ml , M 2 

1,1,0 





M2^ Ml 

1,1,1 






we then defined the ISI-free eapaeity under delay eonstraints, C{d). Afterwards, a suboptimal 
approaeh to determine the eneoder and deeoder funetions (modulation seheme) that gives a lower 
bound on C{d) was proposed. The results show that these simple sehemes are near optimal, i.e, 
beeome elose to the eonstrained eoding eapaeity. 

References 

[1] I. F. Akyildiz, J. M. Jomet, and M. Pierobon, “Nanonetworks: A new frontier in communications,” Communications of 
the ACM, vol. 54, no. 11, pp. 84-89, November 2011. 

[2] T. Nakano, M.J. Moore, F. Wei, A.V. Vasilakos, J. Shuai, “Molecular communication and networking: Opportunities and 
challenges,” IEEE Transactions on NanoBioscience, vol. 11, no. 2, pp. 135-148, 2012. 

[3] Moore, Michael J and Suda, Tatsuya and Oiwa, Kazuhoro, “Molecular communication: modeling noise effects on 
information rate”, IEEE Transactions on NanoBioscience, Vol. 8, No. 2, 2009. 

[4] B. Atakan, S. Galmes, and OB. Akan, “Nanoscale communication with molecular arrays in nanonetworks.” IEEE 
Transactions on NanoBioscience, vol. 11, no. 2, pp. 149-160, 2012. 

[5] M. Pierobon and I. F. Akyildiz, “A statistical-physical model of interference in diffusion-based molecular nanonetworks.” 
IEEE Transaction on Communications, vol. 62, no. 6, 2014. 


DRAFT 


August 12, 2015 




21 


[6] G. Aminian, H. Arjmandi, A. Gohari, M. Nasiri Kenari U. Mitra, “Capacity of Diffusion based Molecular Communication 
Networks in the LTI-Poisson Model,” IEEE International Conference on Communications, June 2015. 

[7] M. S. Kuran, H. B. Yilmaz, T. Tugcu, I. F. Akyildiz, “Modulation techniques for communication via diffusion in 
nanonetworks,” IEEE International Conference on Communications, June 2011. 

[8] A. Einolghozati, M. Sardari, F. Fekri, “Design and Analysis of Wireless Communication Systems Using Diffusion- 
Based Molecular Communication Among Bacteria,”, IEEE Transactions on Wireless Communications, vol.l2, no. 12, 
pp.6096,6105, December 2013. 

[9] H. Arjmandi, A. Gohari, M. Nasiri-Kenari and Farshid Bateni, “Diffusion based nanonetworking: A new modulation 
technique and performance analysis,” IEEE Communications Letters, vol. no. 4, pp. 645 - 648, 2013. 

[10] M. Movahednasab, M. Soleimanifar, A. Gohari, M. Nasiri Kenari U. Mitra, “Adaptive Molecule Transmission Rate for 
Diffusion Based Molecular Communication,” IEEE International Conference on Communications, June 2015. 

[11] Md. Kabir, S. M. Islam, and K.S. Kwak. “D-MoSK Modulation in Molecular Communications.” IEEE Transactions on 
NanoBioscience. 

[12] Pudasaini, Subodh, Seokjoo Shin, and Kyung Sup Kwak. “Run-length aware hybrid modulation scheme for diffusion-based 
molecular communication,” 14th International Symposium on Communications and Information Technologies (ISCIT), 2014. 

[13] Tepekule, Burcu, et al. “Novel Pre-Equalization Method for Molecular Communication via Diffusion in Nanonetworks,” 
IEEE Communications Letters, June 2015. DOI:10.1109/LCOMM.2015.2441726 

[14] M. S. Leeson, D. H. Matthew, “Eorward error correction for molecular communications,” Nano Communication Networks, 
vol. 3, no. 3 , pp. 161-167, 2012. 

[15] P. Y. Ko, et al. “A new paradigm for channel coding in diffusion-hased molecular communications: Molecular coding 
distance function,” Global Communications Conference (GLOBECOM), 2012. 

[16] Shih, Po-Jen, Chia-han Lee, and Ping-Cheng Yeh. “Channel codes for mitigating intersymhol interference in diffusion-based 
molecular communications,” Global Communications Conference (GLOBECOM), 2012. 

[17] B. H. Marcus, R. M. Roth, P. H. Siegel, “An Introduction to Coding for Constrained Systems,” available at https://www. 
math.ubc.ca/~marcus/Handbook 

[18] C. E. Shannon, “A Mathematical Theory of Communication,” The Bell System Technical Journal, 1948. 

[19] G. Bocherer and R. Mathar, “Matching dyadic distrihutions to channels,” IEEE Data Compression Conference (DCC), 

2011 . 


August 12, 2015 


DRAFT 


